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| Abstract. I present a formulation of a Variable Flavour Number Scheme for heavy quarks that is 

implemented up to NNLO in the strong coupling constant and may be used in NNLO global fits for 
f^*) ■ parton distributions. 
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Tj- \ While up, down and strange quarks are treated as effectively massless partons, charm, 

bottom and top have to be regarded as heavy partons. There are two distinct regimes for 
these types of quarks. At low scales, Q 2 ~ m 2 H , they are only created in the final state 
and described using the Fixed Flavour Number Scheme (FFNS) 

£j I F t (x, Q 2 ) = C[f(Q 2 /m 2 H ) ®f k f {Q 2 ). 

However, for Q 2 mjj, the coefficient functions contain large ln(<2 2 / m //) terms, spoil- 
ing the perturbative expansion. In this regime it is more appropriate to treat the quarks 
like massless partons, and the large ln(<2 2 /m^) terms are summed via the DGLAP evo- 
^3 ' lution equations. The simplest recipe involving this regime is the Zero Mass Variable 
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Flavour Number Scheme (ZMVFNS). This ignores all 0(mjf/Q ) corrections, i.e. 

^ F^Q 2 )=Cff F ®f/ + \Q 2 ). 
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The partons in different flavour-number regions are related perturbatively, 

f n k f+ \Q 2 )=A jk (Q 2 /m 2 H )®f'^(Q 2 ), 

where the perturbative matrix elements A jk(Q 2 / m jj) containing ]n(Q /ntjj) terms guar- 
antee the correct evolution for both descriptions. At LO, i.e. zeroth order in as, the 

relationship between the two descriptions is trivial - q(g) k f+ (Q 2 ) = q(g)" k J ' (Q 2 )- At 
NLO, i.e. first order in a s (h+(Q 2 ) = (h + h) (Q 2 )), 

i.e. the heavy flavour evolves from zero at Q 2 = mjj and the gluon loses corresponding 
momentum. It is natural to choose Q 2 = mjj as the transition point. At NNLO, i.e. second 
order in as, there is much more complication 

/f /+ V) = (^) 2 I(4°+4 lln (e 2 /^)+4 2ln2 (e 2 /^))®/; / (e 2 ), 
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FIGURE 1. NNLO F£(x, Q 2 ) and F 2 (x, Q 2 ) in zero-mass VFNS 
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where A -j is generally non-zero |2fl. There is no longer a smooth transition at this order, 
and in fact the heavy parton begins with a negative value at small x. 

This leads to discontinuities in the partons and, without the correct treatment, also in 
the structure functions. ZMVFNS coefficient functions also lead to discontinuities at the 
transition point due to a sudden change in the flavour number in the coefficient functions. 
(This is already true at NLO, i.e. F 2 H (x, Q 2 ) =0 Q 2 <m 2 H ,= ^C 2 , g <g> g n f +l (Q 2 )Q 2 > 
mjj, but the effect is very small.) This is a large effect at NNLO and is also negative at 
smallish x (x ~ 0.001). Hence, ZMVFNS is not really feasible at NNLO, leading to a 
huge discontinuity in Fffa Q 2 ), which is significant in F 2 (x, Q 2 ), as shown in Fig. 1. 

Hence we need a general Variable Flavour Number Scheme (VFNS) interpolating 
between the two well-defined limits of Q 2 < mjj and Q 2 ^> mjj. The VFNS can be 
defined by demanding equivalence of the n y and nf + 1 -flavour descriptions at all orders, 

F^Q 2 ) = Cff(Q 2 /m 2 H )^f n /(Q 2 )=Cfj(Q 2 /m 2 H )^f^ + \Q 2 ) 
= C\j{Q 2 /m 2 H )®A }k {Q 2 /m 2 H )®f n k f {Q 2 ) 
- C[,[(Q 2 /m 2 H ) = ClJ{Q 2 /m 2 H )®A jk {Q 2 /m 2 H ). 

At ^(a s ) this gives 

C F 2 F g \Q 2 /m 2 H ) = C V 2 F H ° H (Q 2 /m 2 H ) ®P>(Q 2 /m|) + C V 2 F g l (Q 2 / m 2 H ) . 

The VFNS coefficient functions tend to the massless limits as Q 2 /mh — > °°, as demon- 
strated to all orders in [3], and if we use the zeroth order cross-section for photon-heavy 
quark scattering we obtain the original ACOT scheme yj]. 



However, ##(<2 2 / m #) * s on ty uniquely defined as Q 2 /m^ — > °°, i.e. one can swap 
^{m 2 H /Q 2 ) tenns between C 2H ^(Q 2 /m^) and C^ ,1 (Q 2 /mj I ). Similar reasoning holds 

for HH(Q 2 / m Jt)- The ACOT prescription violated the threshold W 2 = Q 2 (l -x)/x > 
AM 2 since only one quark was needed in final state. The Thorne-Roberts variable flavour 
number scheme (TR-VFNS) ||4j] recognized this ambiguity and removed it by imposing 
continuity of [dFijd ln<2 2 ) at the transition point. This guaranteed smoothness at Q 2 = 
ntjj, but was complicated and cumbersome when extended to higher orders. 

There have been other alternatives, and most recently the ACOT(^) prescription |0] 
defines F 2 > (x,Q 2 ) = h + (x/x max ,Q 2 ), where x max = Q 2 / '(Q 2 + Am 2 H ) . The coefficient 
functions tend to the massless limit for Q 2 /mj J — > °° but also respect the threshold 
requirement W 2 > Amjj for quark-antiquark production. Moreover it is very simple. For 
the VFNS to remain simple (and physical) at all orders I choose C 2t } H (Q 2 /mj J ,z) = 
C^fjQiz/xmax)- 1 Adopting this convention then at NNLO we have, for example, 

Since A# g (l,z) 7^ 0, C 2 Hg (Q 2 /m H ,z) is discontinuous at <2 2 = m^, and this compen- 
sates exactly for the discontinuity in the heavy flavour parton distribution. 2 
There is one more issue in defining the VFNS: the ordering for F 2 H (x, Q 2 ), i.e. 

n y-flavour nf+1 -flavour 

lo g<g®«»/ cj££®*+ 

NLO (g) (C*g®*»/ + C*g® E»/) g(cjgi®*+ +C 



Switching directly when going from rif to rif + 1 flavours leads to a discontinuity. We 
must decide how to deal with this. Up to now ACOT have used e.g. at NLO 

i.e. the same order of 0(5 above and below, but LO below and NLO above. The Thorne- 
Roberts scheme proposed e.g. at LO 

us(Q 2 ) ^ff,i (Q 2 \^ „, v ^ _ a s (m 2 ) ' ^~ 



-c£*j ( 1 ) »g n '(ml) +Cl r / H » A+ (2 2 



i.e. the higher order «s term is frozen when going upwards through Q 2 = mjj. This 
difference in choice is extremely important at low Q 2 . 



1 It is also important to choose C^^" H (Q 2 jrr^ H ,z) °c C™^(z/x mflx ). 

2 At NNLO there are also contributions due to heavy flavours in loops away from the photon vertex. These 
are included within the VFNS and lead to a discontinuity in the coefficient functions for light flavours 
cancelling that in the light quark distributions. Strictly, part of this contribution should be interpreted as 
light flavour structure functions, while part of it contributes to (x 7 Q 2 ) |8]. 
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FIGURE 2. Comparison of NLO and NNLO predictions for F 2 e [x,Q 2 ) 

Making this choice, in order to define the VFNS at NNLO we need the ^(ttf) heavy 
flavour coefficient functions for Q 2 < mjj. However, these are not yet calculated (making 
a NNLO FFNS problematic). We know the leading threshold logarithms [6], and can 
derive the leading ln(l/x) term from fc^ -dependent impact factors [7], 

C^/° wx {Q 2 /mU) = 96i^/(2 2 />4), /(I) - 4, 

and C 2 F / q l0WX (Q 2 / m HiZ) = 4/9C FF / g lowx (Q 2 /m 2 ,z). By analogy with the known 
NNLO coefficient functions and splitting functions I hypothesize that 

C^ F / g lowx (Q 2 /m 2 H ,z) = y(ln(l/z) -4)(1 -z/x max ) 20 f(Q 2 /m 2 H ), 

i.e. the leading ln( l/z) term is always accompanied by ~ —4, and the effect of the small z 
term is damped as z — > 1. Using the full (if slightly approximate) VFNS one can produce 
NNLO predictions for charm with discontinuous partons, but a continuous F c (x 1 Q 2 ). 
NNLO clearly improves the match to lowest Q 2 data BaLLO], where NLO is generally 
too low, as seen in Fig. 2. 

Hence, I have devised a full NNLO VFNS, with a small amount of necessary mod- 
elling. This seems to improve the fit to the lowest x and Q 2 data greatly. It also guarantees 
continuity of the physical observables, such as structure functions, despite the disconti- 
nuity in NNLO parton distributions. It can now be used in a full NNLO global analysis. 
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